Multiple Thresholds in a Model System of Noisy Ion Channels 
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Voltage-activated ion channels vary randomly between open and closed states, influenced by the 
membrane potential and other factors. Signal transduction is enhanced by noise in a simple ion 
channel model. The enhancement occurs in a finite range of signals; the range can be extended 
using populations of channels. The range increases more rapidly in multiple-threshold channel 
populations than in single-threshold populations. The diversity of ion channels may thus be present 
as a strategy to reduce the metabolic costs of handling a broad class of electrochemical signals. 
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Voltage-activated ion channels are essential elements 
in biological signal transduction, playing important roles 
in synaptic transmission, generation of neural action po- 
tentials, regulation of membrane potentials and intracel- 
lular Ca^"*" concentrations, and other cellular functions 
P, 0, 0, H . The gating dynamics of the channels allow 
the nonconductive cell membrane to take part in electri- 
cal conduction and signaling. Channels vary between a 
conducting or open state and a nonconducting or closed 
state, with intermediary states in the transition being un- 
stable and short-lived. The transition between open and 
closed states is influenced by a broad assortment of fac- 
tors, principally the membrane potential, but also includ- 
ing hormones, toxins, protein kinases and phosphatases, 
and thermal fluctuations. Volt age- activated channels are 
functionally diverse in their sensitivity to depolarization; 
indeed, Lee et al. identify no fewer than five distinct 
activation thresholds for Ca^"*" channels. 

Channel gating dynamics are intrinsic stochastic tran- 
sitions that depend strongly on external factors, so that 
the channel constitutes a single-molecule sensor or com- 
munication channel, transforming membrane potentials 
into noisy ionic currents. Noise can have surprising ef- 
fects in many nonlinear systems. Perhaps the best known 
of these is stochastic resonance (SR), wherein the pres- 
ence of noise enhances the response of a thresholding sys- 
tem to a weak periodic signal (for a review, see Ref. 
SR has been experimentally demonstrated in a system of 
parallel ion channels |^ , and studied theoretically in nu- 
merous systems (see, e.g., Refs. ^J3„,1Q3 JT, 12, 13, 14]). 
In this work, we use a simple ion channel model to fur- 
ther investigate the stochastic nature of voltage- activated 
ion channels, characterizing the limits that noise imposes 
on information transduction in individual channels and 
systems of channels. 
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We make use of a discrete model in which the channel 
switches between distinct open and closed states, omit- 
ting the dynamics of the transition process. Such a dis- 
crete model captures the bistable nature of the channel 
dynamics. By omitting the transition dynamics, we as- 
sume that any stimulus to the channel varies slowly com- 
pared to the time scale of channel opening and closing. 
We can extend this assumption to a quasistatic approxi- 
mation, where the channel is always in equilibrium, and 
describe the channel opening probability by using the 
steady state (time t — > oo) probability for the permis- 
sive state also called the activation function of the 
channel. This probability is given by 
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where z is the (dimensionless) valence of the "gating par- 
ticles," F is Faraday's constant, V is the transmembrane 
potential, R is the ideal gas constant, and T is the tem- 
perature. The parameter Vb is a bias (or noisy threshold) 
in the potential to which the channel tends to open. Gal- 
ley et al. [13 have investigated noise in model ion chan- 
nels essentially identical to the ones defined by eq. 
demonstrating that stochastic resonance occurs. 

For notational clarity, we lump several of the parame- 
ters into a thermal noise parameter a, such that 
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The definition in eq. Q amounts to changing the units 
of measure for the temperature to Volts. At room tem- 
perature and with 2: = 1, a is approximately 25 mV. 
Using eq. (0), the channel opening probability becomes 



p(F;a) = 
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with a corresponding probability of remaining closed of 
q{V]a) — 1 —p(V]a). For single channels, or popu- 
lations of channels with homogeneous behavior, we can 
take Vo — without loss of generality; the behavior for 
other values of Vq can be recovered by translating the 



2 



potential in eq. © by the desired value for the thresh- 
old potential. In heterogeneous populations of channels, 
the differing values of the threshold potentials can have a 
profound effect on the behavior of the system of channels, 
as we will demonstrate below. 

For a system of N channels, each of which are exposed 
to the same transmembrane potential V but with differ- 
ent realizations of the noise (i.e., we have i.i.d. noise), we 
let Zn be the number of channels that open during the 
time interval. The expected state of the membrane, i.e., 
Zf^, can be calculated in a straightforward fashion. The 
expectation value {Zj\f) (V; a) and variance {V; a) 
can be expressed using p {V; a), q {V; a), and N, giving 

(Zjv)(F;a) = Np{V;a) (4) 
a|„(F;a) = Np {V; a) q {V; a) . (5) 

In calculating eqs. Q and we have made use of the 
independence of the noise for the channels. 

We will focus on the behavior of ion channels near 
the threshold value. To explore the ability of the ion 
channels to serve as a transducer of electrical signals, 
we reproduce the input potential by decoding the state 
of the membrane (i.e., the numbers of open and closed 
channels). Near the threshold, this gives rise to linear 
decoding rules. The basic approach is similar to the "re- 
verse reconstruction" using linear filtering that has been 
applied with great effect to the ana.lysis of a number of 
biological systems (see, e.g., Refs. [S 113 111 III Ei). 
Despite the quite different nature of the signals we con- 
sider, our static population analysis (below) shares the 
key features of the temporal analyses in the cited works. 
Specifically, by comparing the actual input to an esti- 
mate derived from the systemic response, we can work at 
greater remove from the details of the systems. We thus 
bypass the need for complete descriptions of the trans- 
mission process and of the manner in which the system 
response is used. This permits general limits to be estab- 
lished on the signal-transducing ability of a population of 
channels, regardless of what the actual mechanisms may 
be. 

We expand the expected number of open channels 
(Zpf) (V; a) to first order near the threshold (i.e., V ^ 0), 
giving 



{Zn) {V;a) = ^ + ^V + 0{V^) 
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An example of the linear approximation is shown in fig.^ 
Dropping the higher order terms and inverting eq. © 
suggests a linear decoding rule of the form 

^„.4»(f -i) . ,7, 

where Vn is the estimate of the input potential. Com- 
bining eqs. Q and Q, we can show that 



VN){V;a)^4a[p{V;a)-^] . (8) 
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FIG. 1: First order approximation of the expectation value for 
the channel opening. Near the threshold potential (Vb = 0) 
of the channel, the expectation is nearly linear. Further from 
the threshold, the expectation value saturates at either zero 
or one, where the linear approximation diverges from the true 
value. The values shown here are based on thermal noise 
a = 1, giving a response width W = 2. 



The expected value of Vn is thus seen to be independent 
of N; for notational simplicity, we drop the subscript and 
write {V). Note that, as the thermal noise increases, the 
expected value of the decoded potential closely matches 
the input potential over a broader range. 

We must also consider the uncertainty of the potential 
value decoded from the state of the membrane. This 
leads to a total decoding error AV^v with the form 
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Using eq. 0, we can derive several properties that are 
useful for understanding the role of noise in the channel 
behavior. In particular, it is straightforward to show that 
{V) {V; a) = V{V) (1; {a/V)), e {V; a) = Ve (1; {a/V)), 



and CT? (V; a) 

Vn 



V^al {l;a/V), for all V 0. Thus, 

Vjv 



the noise dependence for both the reconstructed stimu- 
lus and its uncertainty can be understood with a single 
stimulus. The total error AV^ {1; (a/V)) is minimized 
for a nonzero value of the noise parameter a, analogous 
to the stochastic resonance effect; see fig. |21 

In fig. 121 we show how AV^ varies with the number of 
channels N. As N increases, the error curve flattens out 
into a broad range of similar values. Thus, the presence 
of noise enhances signal transduction without requiring 
a precise relation between V and a. This effect is anal- 
ogous to the "stochastic resonance without tuning" first 
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FIG. 2: Comparison of decoding error sources. The values 
shown here are calculated for a single channel. The mini- 
mum in occurs for nonzero thermal noise, analogous to 
stochastic resonance. 
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FIG. 3: Effect of the number of channels on the decoding 
error. As A'' becomes large, the error curve flattens out, in- 
dicating a broad range of noise values that all give similar 
accuracy in the decoding process. 

reported by Collins et al. '2l1 and previously explored in 
two-state systems not dissimilar to the one investigated 
here (see, e.g., Ref. H^)- 

Informally stated, SR without tuning allows for a wider 
range of potentials to be accurately decoded from the 
channel states for any particular value of a. To make 
this notion of "wider range" precise, we again focus our 
attention on the expected behavior of the channels (see 
fig-D- The expected channel response {Z) matches well 
with the linear approximation when \y\ < a. From this, 
the width W can be defined to be 2a. Other definitions 
for the response width are, of course, possible, but we still 
should observe that the width is proportional to a, since 
the probability for channel opening depends only on the 
ratio of V and a (eq. Q). The same width is found for 
multiple identical channels, because the total expected 
current is proportional to the single channel behavior, 
without broadening the curve in fig. ^ 

The response width can thus be increased by increas- 
ing the thermal noise parameter a. As seen in figs. |5| 
and 13 such an increase ultimately leads to a growth in 



the decoding error AV)^. As a becomes large, AV)^ is 
dominated by cr? and we have the asymptotic behavior 

VjV 

^Vl{y-a) = o(^ . (12) 

based on eq. Hll|) . The growth in AV^ with increas- 
ing a thus can be overcome by further increasing the 
number of channels in the population. Therefore, the re- 
sponse width W is effectively constrained by the number 
of channels N , with W = 0{^/N) for large N. 

An arbitrary response width can be produced by as- 
sembling enough channels. However, this approach is 
inefficient, and greater width increases can be achieved 
with the same number of channels. Consider instead di- 
viding up the total width into M subranges. These sub- 
ranges can each be covered by a subpopulation of N chan- 
nels, with the subpopulations having different thresh- 
olds from one another. The width of each subrange is 
0{VN), but the total width is 0{MVN). Thus, the 
total response width can increase more rapidly as addi- 
tional types of channels are added. Conceptually, multi- 
ple types of channels arise naturally as a way to provide 
a wide range of accurate responses, with multiple chan- 
nels in each type providing independence from any need 
to "tune" the thermal noise to a particular level. 

To describe the behavior of channels with different 
thresholds, much of the preceding analysis can be directly 
applied by translating the functions along the potential 
axis to obtain the desired threshold. However, system be- 
havior was previously explored near the threshold value, 
but heterogeneous populations of channels have multi- 
ple thresholds. Nonetheless, we can produce a compara- 
ble system by simply assessing system behavior near the 
center of the total response width. 

To enable a clean comparison, we set the thresholds in 
the heterogeneous populations so that a linear decoding 
rule can be produced. A simple approach that achieves 
this is to space the thresholds of the subpopulations by 
2W = 4a, with all channels being otherwise equal. The 
subpopulations with lower thresholds provide an upward 
shift in the expected number of open channels for higher 
threshold subpopulations, such that the different subpop- 
ulations are all approximated to first order by the same 
line. Thus, the expected total number of open channels 
leads to a linear decoding rule by expanding to first or- 
der and inverting, as was done earlier for homogeneous 
populations. Note that this construction requires no ad- 
ditional assumptions about how the channel states are to 
be interpreted. 

To illustrate the effect of multiple types of channels, we 
begin with a homogeneous baseline population M — 1 oi 
N — 1000 channels with Vq = and apply a potential 
V with thermal noise a = 1. Using the definition above, 
the response width is 14^ = 4. We then consider two 
cases, homogeneous and heterogeneous, in each of which 
we increase the response width by doubling the number 
of channels while maintaining similar error expectations 
for the decoded currents. 
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FIG. 4: Total decoding error in homogeneous and heteroge- 
neous populations of channels. The heterogeneous channel 
population (M — 2,N = 1000) has a broader "basin" of low 
error values than the baseline (M = 1, A'' = 1000) and homo- 
geneous (M = 1, A'^ = 2000) populations. 



In the homogeneous case, we have a single population 
(Af — 1) with N — 2000 channels. Doubling the num- 
ber of channels allows us to increase the temperature 
parameter a by a factor of a/2 with similar expected er- 
rors outside the response width. Thus, we observe an 
extended range, relative to the baseline case, in which 
we can reconstruct the stimulus potential from the state 
of the channels (fig. EI . 

In the heterogeneous case, we instead construct two 
subpopulations (M = 2) with N ~ 1000 channels. We 
leave the thermal noise parameter unchanged at a = 1. 
One of the subpopulations is modified so that the chan- 
nel thresholds lie at +W ~ 2, while the other is modified 
so that the channel thresholds lie at —W = —2. The 
resulting system of channels has a broad range in which 
we can reconstruct the stimulus potential with low er- 
ror, markedly superior to the baseline and homogeneous 
cases (fig.EJ. The approach used in this example can be 
directly extended to three or more subpopulations. 

The foregoing analysis and example suggest that the 
diversity of channel types found in a living cell are present 
as an information processing strategy, providing a means 
to efficiently handle a broad class of electrochemical sig- 
nals. The superior scaling properties of heterogeneous 
populations of channels can have a profound impact on 
the cellular metabolism; large numbers of channels imply 
a large energetic investment, both in terms of the proteins 
needed to construct the channels and of the increased de- 
mand on ion pumps that accompanies the greater ionic 



currents |3. The action potentials generated in neu- 
rons can require a significant energetic cost [S^], making 
the tradeoff between reliably coding information and the 
metabolic costs potentially quite important. 

In this picture, we expect that different types of cells 
will require different numbers of functionally different ion 
channels. Cells that perform sophisticated signaling and 
respond to a broad variety of signals will need a large 
number of functionally different ion channels, while cells 
that are more specialized to a narrower class of signals are 
likely to have a smaller number of functional types. This 
appears to be generally consistent with the comparatively 
large varietyof ion channels found in excitable cells such 
as neurons 

It is interesting to reconsider the work of Lee et al. 
in light of the analysis presented here. In their work, 
five Ca^"*" channels with distinct thresholds are tabulated, 
including three T-type channels interpreted as being im- 
portant for the electrical responsiveness of neurons. How- 
ever, a full understanding of the physiological roles of 
the channels remains to be found. The present analysis 
indicates that signal transduction by any given type of 
channel is intrinsically limited, regardless of the details 
of how the state of the channels is used. The multiple 
types of channels with their various thresholds provide a 
metabolically and evolutionarily favorable means to over- 
come those limits. 

Although we have used a specific model consisting of 
channels with thermal fluctuations modulating an input 
potential, we expect that the key result is more widely 
applicable. The demonstration of the advantage of multi- 
ple channel types largely follows from two factors that are 
not specific to the model channels. First, the distance of 
the input potential from the threshold is proportional to 
the level of the thermal noise, and, second, the total vari- 
ance of the inputs to the channels is proportional to the 
number of channels. Ultimately, a multiplicity of func- 
tional types of channels with varying thresholds arises 
because the independently distributed noise provides a 
natural scale for the system. 
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